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We employed ab-initio molecular dynamics to directly simulate the effects of Ag alloying (~ 5% Ag con- 
centration) on the phase change properties of Ge2Sb2Te5. The short range order is preserved, whereas a 
slight improvement in the chemical order is observed. A slight decrease in the fraction of tetrahedral Ge 
(sp 3 bonding) is reflected in the reduction of the optical band gap and in the increased dielectric constant. 
Simulations of the amorphous to crystalline phase change cycle revealed the fact that the crystallization 
speed in Ago.sGe2Sb2Te5 is no less than that in Ge2Sb2Tes. Moreover, the smaller density difference and 
the larger energy difference between the two phases of Ago.5Ge2Sb2Te5 (compared to Ge2Sb2Tes) suggest a 
smaller residual stress in devices due to phase transition and improved thermal stability for Ago.sGe2Sb2Te5. 
The potential viability of this material suggests the need for a wide exploration of alternative phase change 
memory materials. 



The computational design of materials is still in its 
nascent stages, but is widely recognized to be one of 
the prime frontiers of materials science. The challenges 
are daunting for several reasons, among these: time and 
length scales drastically different in simulation compared 
to laboratory samples, the need for realistic interatomic 
interactions (nowadays largely based upon pseudopoten- 
tials and density functional theory) leads to tremendous 
demand for computational resources. In the case of the 
phase change memory materials, with compositions near 
Ge2Sb2Te5 (GST), there is clear evidence that current 
first principles simulations can accurately simulate phase 
changes on the time scales accessible to these codes 1 . 
Other work suggests that key quantities like crystal- 
lization speed can be meaningfully inferred from such 
simulations^. These materials are of great fundamental 
interest for their ultrafast phase changes and are the lead- 
ing candidate to replace current non-volatile computer 
memory, a multi-billion dollar market. 

In this paper, we explore new candidate phase change 
materials and show that a silver-doped variant may be 
superior to conventional GST. We elucidate the process 
of crystallization in atomistic detail and particularly note 
the role of the Ag in producing more stable and chem- 
ically ordered materials. Beside the specific prediction 
that the Ag alloy systems deserve careful exploration, we 
highlight the existence of a promising unexplored which 
strongly suggests that others probably exist as well, and 
deserve exploration. 

We have implemented an ab-initio molecular dy- 
namic (AIMD) simulations to study the ultrafast crys- 
tallization of Ag-doped (alloyed) Ge2Sb2Te5. The AIMD 
calculations were performed using Vienna Ab-initio 
Simulation Package (VASP)2r— to generate models of 
Ge 2 Sb 2 Te 5 and Ago.sGesSbsTes (AGST) with 108(24 
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Ge atoms, 24 Sb atoms and 60 Te atoms) and 114(24 
Ge atoms, 24 Sb atoms, 60 Te atoms, and 6 Ag atoms) 
atoms, respectively. The calculations were performed by 
using the projector augmented- wave (PAW) 5,6 method to 
describe electron-ion interactions. The Perdew-Burke- 
Ernzerhof (PBE}£ exchange correlation functional was 
used throughout. Molecular-dynamics (MD) simulations 
were performed in a cubic supercell with a time step of 
5.0 fs using periodic boundary conditions at constant vol- 
ume for annealing, equilibrating and cooling, whereas, 
zero pressure conjugate gradient (CG) simulations were 
performed for relaxation. The final models were prepared 
by using the "Melt and Quench" method 8 starting with a 
random configuration at 3000K. Densities of 6.046 gem -3 
and 6.234 gem -3 , respectively for GST and Ag-GST, 
were used during the process. After mixing the random 
configurations at 3000K for 20ps, each model was cooled 
to 1200K in lOps and equilibrated for 60ps. A cooling 
rate of 12K/ps was adopted to obtain the amorphous 
models from the melt at 1200K to 300K and followed by 
equilibration at 300K for another 50ps. Finally the sys- 
tems were fully relaxed to a local minimum at pressure. 
Each of these models was then equilibrated at 300K for 
25ps and data was accumulated for the last lOps and 
statistically averaged to study the structural properties. 
Three independent models were generated for both struc- 
tures. Unless otherwise stated the results presented are 
the statistical average of the models. 

The zero-pressure structures of amorphous and crys- 
talline GST and AGST are presented in FigHJ Inter- 
estingly, the computation of the atomic densities shows 
a relatively small density change (4.61%) between the 
two phases of AGST in contrast to a density change 
of 6.84 % in pure GST. This smaller volume (density) 
change in Ag-doped GST could result in reduced resid- 
ual stress in PCM devices. We further computed the dif- 
ference in the energies between the amorphous and the 
crystalline phases in GST and AGST. The energy differ- 
ence of 80 meV/atom in AGST is about 20meV/atom 
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FIG. 1. (Color online) Simulated structures (relaxed) of a) 
a-Ge 2 Sb 2 Te 5 , b) a-Ago.5Ge 2 Sb 2 Te5 , c)c-Ge 2 Sb 2 Te 5 , and 
d) c-Ag . 5 Ge 2 Sb 2 Te 5 . (Model 3). Color code; Orange- Te, 
Green- Ge, Purple- Sb, and Gray-Ag. 
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FIG. 2. (Color online) Ge-Te partial pair correlation func- 
tions (PPCFs) at 300K in both phases of Ge 2 Sb 2 Te 5 and 
Ago.5Ge 2 Sb 2 Te5. Only one model (Model 3) from each of the 
Ge 2 Sb 2 Te5 and Ago.5Ge 2 Sb 2 Te5 models is presented for the 
illustrations purpose. 



more than that of pure GST. This larger energy differ- 
ence might yield better thermal stability in Ag-doped 
GST, and could improve the data retention capability 
of PCM devices. To investigate the Ag-induced modifi- 
cation of GST network, we analyzed the local structure 
around Ge via partial pair correlation functions (PPCFs) 
and bond angle distributions (BAD). We inspected the 
Ge-Te PPCF (figEJ) and Ge-centered BAD (figEj) because 
Ge atoms undergo a dramatic change upon phase tran- 
sition i.e. perfect octahedral configurations (p-bonding) 
in the crystalline phase and tetrahedral geometry(sp 3 - 
bonding) in the amorphous phased. In the amorphous 
phase, the Ge-Te bond length is found to be increased 
(by 0.02A) in Ag-doped GST. This Ag-induced change is 
also observed in the Ge-centered BAD as a suppression 
near 109° depicting a reduction in the fraction of tetra- 
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FIG. 3. (Color online) Ge-centered Bond angle distribu- 
tions (BADs) (Model 3) in both phases of Ge 2 Sb 2 Tes and 
Ago. 5 Ge 2 Sb 2 Te 5 at 300K. 



hedral Ge sites due to addition of Ag. This reduction 
caused the increase in the average Ge-Te bond length 
since Ge-Te bond lengths with tetrahedral geometry is 
smaller than Ge-Te with octahedral geometr y 1Q i n . In 
contrast to the amorphous phase, we observed identical 
Ge-Te bond lengths between GST and AGST in the crys- 
talline phases. This is consistent with the fact that all 
the tetrahedral Ge changes to an octahedral geometry 
during crystallization. On the other hand, Ag-induced 
change is negligible in the Sb-Te PPCF and Sb-centered 
BAD, in both phases. This is also supported by the fact 
that Sb always takes octahedral geometry in either of 
the phases. Beside Ge-Te and Sb-Te bonds pairs, we ob- 
served a significant fraction of Ge-Ge, Ge-Sb, Sb-Sb and 
Te-Te bond pairs as nearest neighbors especially in amor- 
phous phase. These "wrong bonds"—, amount to 25% in 
amorphous phases and fall to about 6% in the crystalline 
phases. From analysis of the local structure, we are also 
able to identify the interaction of the dopants (Ag atoms) 
in the host network. The Ag PPCF confirms that Ag is 
mainly bonded to Te rather than to Ge or Sb. This is 
also true in crystalline phase where Ag takes the vacancy 
sites (or similar sites as Ge/Sb). 

The investigation of the electrical properties [via elec- 
tronic density of states(EDOS)] in both phases of GST 
and AGST confirms no major differences in the EDOS, 
with p-like states of Te, Sb and Ge dominating both the 
valence and the conduction band and Ag contributing a 
d-like state about 4eV below the fermi level. The band 
gap is observed to decrease with the presence of Ag. Since 
the larger band gap in a-GST as compared to c-GST is 
due to the presence of sp 3 -bonded Ge atomsi^, the re- 
duced band gap by doping can also be attributed to the 
reduction of the tetrahedral Ge atoms. 

The utility of the PCMs stems from the contrast in 
optical properties between amorphous and crystalline 
phases. The imaginary part and the real part of the di- 
electric function confirm that the optical contrast is pre- 
served in AGST. These results are in consistent with the 
results reported by Shportko et al*&. The estimation of 
the optical dielectric constant i.e. the lower energy-limit 
of the real part of the dielectric function (uo — )• 0) is pre- 
sented in Table H We observed a slightly higher dielectric 
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TABLE I. Comparison of dielectric constant between the two 
phases of Ge 2 Sb 2 Te 5 and Ago.5Ge 2 Sb 2 Te5. (Model 3) 



Material 


Amorphous 


Crystalline 


%increase 


Ge 2 Sb 2 Te 5 


25.9 


53.0 


105 


Ge 2 Sb 2 Te 5 (Ref^) 


16.0 


33.3 


108 


Ago. 5 Ge 2 Sb 2 Te 5 


26.9 


60.2 


124 



constant in AGST as compared to GST and suspect that 
this is due to improved medium-range order (increase 
in the number of four-membered, near-square, rings^) in 
AGST. 

The full potential of simulation is revealed in directly 
simulating phase transitions^^. We annealed the a-GST 
and a-AGST models at 650K until each of the models 
crystallized. The process proceeds in three steps (I, II, 
and III), as explained by Lee et al 15 . Period I is termed 
the incubation period. Period II is the main time seg- 
ment in which the process of crystallization occurs and 
the third Period (III) defines the completely crystallized 
state. To understand the crystallization process we ob- 
served the evolution of total energy of the system, the 
number of 4- member rings, seeds in the spirit of clas- 
sical nucleation theory (CNT) 15 , the number of wrong 
bonds, and the total coordination numbers as a func- 
tion of time and present in FigHJa-d). We observe al- 
most no change in the total energy and the number of 
4- member rings during the incubation period (I), how- 
ever, we observe a significant decrease in the number of 
wrong bonds. Wrong bonds keep declining during the 
crystallization period (II) until the crystallization is com- 
plete. About 5 to 7 % of wrong bonds(mainly Ge-Sb 
and Sb-Sb bonds) still exist even after the crystalliza- 
tion is occurred. The total energy and the number of 
4-member rings on the other hand are found to be cor- 
related to each other, with the number of rings increas- 
ing monotonously during the crystallization period. We 
further computed the evolution of pair correlation func- 
tions and the Ge-centered bond angle distribution (for 
model 3) and present these findings in figures [5] and [6] 
respectively. The top panels of FigJS] represent the to- 
tal pair correlation functions (TPCF). The middle panel 
on the left represents the X-Te (where X=Ge and Sb) 
pair correlation functions whereas that on the right rep- 
resents X-Te (X=Ge, Sb, and Ag) PCFs. Finally, the 
bottom panels represent the correlation of wrong bonds 
(Y-Y, Y = Ge, Sb, or Ag). These figures clearly depict 
the evolution of medium to long range order (secondary 
peaks in PCFs) which is the signature of the crystalline 
structures. The prominent medium to long range order 
peaks start evolving during the Period II. The average 
peak positions at 2. 98 A, 5. 3 A, and 6. 8 A for X-Te and 
4.2A and 7.4A for wrong bonds well represent the crys- 
talline GST structure. Similarly, the Ge-centered BAD 
(FigEJ) shows an evolution of narrow and prominent dis- 
tribution around 90° and 180° during the Period II. The 
narrowing of the peak at 90° illustrates the conversion 
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FIG. 4. (Color online) Comparison of the total energy (a), 
the number of four- fold rings (b), the number of wrong bonds 
(c), and the total coordination numbers (d) as functions of 
time in Ge2Sb2Tes (dark, black) and Ago.5Ge2Sb2Te5(light, 
green). (Model 3). The vertical lines separate the three peri- 
ods. 



of the tetrahedral Ge (angular distribution at 109°) into 
the octahedral Ge. The peak at around 180° becomes 
visible during the period II where the total coordination 
numbers (CN) reaches about 4.5, similar to the evolution 
of secondary peaks in PPCFs. CN on the other hand 
depicts a correlation with the total energy of the system 
i.e. CN is almost constant during the incubation period 
(I), increases during the crystallization period (II), and 
becomes constant after the crystallization is established. 

To investigate the effect of Ag on the crystallization we 
compare total energy, the number of four member rings 
and the coordination number (FigHJ. Since the crys- 
tallization of three different models of pure GST shows 
large fluctuations in the duration of Period I and II, es- 
pecially Period I, the estimation of crystallization time 
shows significant uncertainty. The incubation periods 
(Period I) in three different pure GST models vary from 
50 ps to 200ps whereas the crystallization periods (Pe- 
riod II) vary from 40ps to 150ps. These times in AGST 
are (80-1 lOps) for incubation periods and (70-1 lOps) for 
crystallization periods. To understand this we examined 
the local structures of the starting configuration of the 
three GST models. Unsurprisingly, we observe a clear 
distinction in the number of wrong bonds and four mem- 
bered rings, i.e. the model with the short incubation and 
crystallization period has the least number of wrong bond 
and most number of four membered rings (more ordered 
in a plane). To compare the crystallization of GST and 
AGST we chose configurations with similar numbers of 
the wrong bonds. We observed a clear contrast in the 
duration of both the incubation period and the crystal- 
lization period in these two networks. Both of the periods 
were shorter in AGST than in GST. Total of these two pe- 
riods in AGST measure about 200ps against about 315ps 





suggests that Ago.5Ge2Sb2Te5 is thermally more stable 
than Ge2Sb2Te5. On the other hand, smaller density dif- 
ference in Ago.5Ge2Sb2Te5 between the two phases as 
compared to Ge2Sb2Te5 could well reduced the resid- 
ual stress in the PCM devices. Furthermore, the in- 
creased optical contrast between the two phases as well 
as a potential increase in crystallization speed might lead 
to PCM devices with improved performance. 
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FIG. 6. (Color online) Time evolution of Ge-centered 
bond angle distributions in Ge2Sb2Tes (left) and 
Ago.sGesSbsTes (right). (Model 3). 



in pure GST, clearly suggesting a faster crystallization in 
Ag-doped GST. 

In conclusion, we have used AIMD simulations to study 
the effect of Ag doping in Ge2Sb2Tes and directly sim- 
ulate ultrafast phase transitions. Medium range order 
is found to be improved with the addition of Ag in 
the form of increased number of four membered rings 
and decreased fraction of tetrahedral Ge. We were also 
able to simulate the process of amorphous to crystalline 
phase transition. The incubation and crystallization pe- 
riod were found to depend on the wrong bonds present 
in the amorphous phase. Moreover, our simulation re- 
vealed that the crystallization speed is increased by dop- 
ing Ge2Sb2Te 5 with Ag. The larger energy /atom dif- 
ference between amorphous and crystalline phases also 
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